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Abstract 

The stretching of brushes of long polymers grafted to a planar surface is investigated by Monte Carlo sim- 
ulations in the limit of very high grafting densities, as achieved in recent experiments. The monomer density 
profiles are shown to deviate considerably from the paraboUc Umiting form predicted by self-consistent field 
theory. A rapid transition is observed from parabolic to fully stretched polymers, characterized by a dra- 
matic change in the end-monomer height distribution and by a clear cross-over in the slope of the brush 
height versus scaled grafting density. 

PACS numbers: 68.47 .Mn, 68.47.Pe, 61.25.Hq, 82.35.Lr 
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There have been numerous experimental [[D, theoretical [l2l[3l|4l and simulation [151[6l studies 
of polymer "brushes" grafted to substrates, confirming the Alexander scaling of the brush height 
with grafting density [7J. They show that the monomer density profiles go over to a universal 
parabolic profile as a function of the distance from the surface, for sufficiently strong stretch- 
ing, as predicted by self-consistent field theory (SCF) [|3l|4l[8l. Recent experimental techniques 
based on surface-initiated polymerization BH [TOl allow higher grafting densities than previously 
achieved. In this letter we analyse the break-down of SCF theory for long chains and high grafting 
densities by advanced Monte Carlo simulations of polymers in good solvent. Upon increasing 
the grafting density, we observe a transition from the parabolic regime to fully stretched brushes 
which correlates with a dramatic change in the end-monomer density profile. The exposure of the 
end monomers of very dense brushes opens new possibilities for the development of chemically 
active soft surfaces. 

Consider a layer of N identical flexible polymers of M monomers and contour length L = 
(M — l)b (where b is the segment length), in good solvent. The first monomer of each chain is 
chemically grafted to a grafting site on a planar substrate of area A; the dimensionless grafting 
density is defined as a = Nb^/A. For convenience b will henceforth be the unit length. The key 
length scales are the mean distance between grafting sites, d ~ o^^/^, the radius of gyration of an 
isolated polymer, Rg ~ V (where the Flory exponent v = 0.588 under good solvent conditions) 
and the brush height h which, according to the Alexander-de Gennes mean field argument [l2l|71, 
scales like 

h Lo^ ~ Lo°-3^ ~ Ld-^-^ (1) 

i.e. the chains stretch upon increasing the grafting density due to their mutual repulsion. This 
scaling has been confirmed by neutron scattering diffraction experiments yj. An important di- 
mensionless parameter is Og = oR^ R^/d^. When < 1, polymers grafted to neighbouring 
sites do not overlap and the height h of quasi-independent coils scales like /? ~ i?g ~ ("mush- 
room regime"). When Og > I, neighbouring coils interact and begin to stretch, i.e. h L. This 
regime has been widely studied within SCF theory [|2l [3 [8l and by MC simulations of various 
models, involving in general rather short chains (L < 100) [|51|6l[8l[IIl[l2l. The polymer "brush" 
may be characterized by the monomer distribution function P{z), where z is the vertical distance 
from the substrate, normalized such that Jq P{z)dz = N. Within SCF, the dimensionless control 
parameter is (3 = {h/RgoY (where Rgo ~ L^^^ is the radius of gyration of the non-interacting, Gaus- 
sian chain), and in the strong stretching regime (3 >> 1, SCF reduces to its "classical" limit, which 
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Figure 1: (Color on-line)Scaled monomer distribution function versus scaled altitude for brushes of SAW 
polymers of length L = 50, and several grafting densities 0.002 < a < 0.5. The two lowest a correspond to 
the "mushroom" regime, while deviations from SCF scaling occur for a > 0.25. The inset shows the end 
monomer altitude profiles versus z/L for the same values of a 

predicts a universal parabolic curve when P{z)o^^^^^^^ is plotted versus z/ (^Lo^^^^^^^^)j [[31 IH. 
This scaling is confirmed by MC simulations, with deviations from the parabolic shape at small 
and large z in the case of moderate stretching (H [6l [8l [HI [I2l. Some deviations from a universal 
master curve, due to a flattening of the profiles at intermediate heights, have been observed in sim- 
ulations at the highest grafting densities Og explored so far [[8l[T3l[l4l. Meanwhile experimental 
grafting techniques based on surface-initiated polymerization BUl have achieved extremely high 
grafting densities, up to nearly one site per nnP' [|9l[T0ll. 

To explore the structure of polymer brushes in the regime of such high grafting densities, cor- 
responding to the regime of Og >> 1, which may be achieved by increasing the grafting density o 
or considering very long polymers, we have carried out extensive MC simulations of the simplest 
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lattice model of flexible polymers in good solvent, namely self (and mutually) avoiding walks 
(SAW) on a cubic lattice. The grafting sites are on a square lattice of identical spacing b. We used 
"annealed" grafting conditions whereby the first monomer of each chain can jump between near- 
est neighbor grafting sites in order to speed up equilibration [fT2l. Three types of MC moves were 
used: discrete translations of entire chains on the square lattice in the x — y grafting plane; chain 
re-growth moves based on the Configurational Bias Monte Carlo (CBMC) algorithm ifTSl ; seg- 
ment re-growth moves based on CBMC with fixed end-points [fTSll . This combination of moves is 
expected to guarantee good equilibration of the simulated brush up to very high grafting densities 
(corresponding to moderately high monomer volume fractions in the brush). The initial conditions 
were generated with fully stretched polymer conformations attached to randomly chosen grafting 
centres on the square lattice. The grafting centres were allowed to diffuse on the square lattice 
(annealed samphng), which is equivalent to a random sampling of fixed grafting points, followed 
by a statistical average over the randomly chosen initial configuration [fT2ll . After initial equili- 
bration, statistical averages of the monomer and end-monomer profiles were taken over typically 
2-10^ MC moves. 

We have systematically computed the monomer and end-monomer distribution functions P{z) 
and Pm{z) over a wide range of a, and for lengths L =50, 200, 400 and 800. Scaled profiles for 
short chains (L = 50) and several o are shown in Fig.[T| The universal scaling regime is reached 
for o > 0.04 (Og > 0.6), and persists up to o = 0.25 (Og ~ 4). Below o ^ 0.04 (Og < 0.6), 
interactions between neighbouring grafted polymers become negligible so that the "mushroom" 
regime is reached, and the profiles no longer follow the SCF scaling. The end-monomer profile 
Pm{z) moves to higher z as o increases, as one might expect. The probabihty of the end monomer 
"returning" towards the grafting surface decreases with increasing o [[8]|. Fig. [2] shows similar 
results, but for longer chains (L =200). The mushroom regime is only observed for o < 0.004 
(Og < 0.3), while the scaling regime is satisfied for 0.01 < a < 0.05 (0.8 < Og < 4). For a > 
0.05, the profiles deviate increasingly from the scaling prediction, and for o > 0.2(Og > 16) P{z) 
flattens and stretches considerably. Concomitantly the end- monomer profiles Pm{z) move further 
and further towards their upper limit, z = L. These trends are considerably enhanced in the case 
L = 800, as shown in Fig. |3j The SCF scaling regime only holds for o < 0.02 (Og < 10) while 
for higher a, P{z) flattens and stretches towards the rectangular profile (P{z) = 1, z < L; P{z) = 0, 
z> L) which is the exact limit for the lattice model, when o = 1 (each polymer then reduces to 
a rod of length L because of the excluded volume constraint). This strong stretching behavior 
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correlates well with the behavior of Pm{z), which undergoes a dramatic transition around a = 0.1 
(Og = 40), when the peak in the end-monomer distribution function sharpens and moves towards 
z = L (in the limit a ^ 1, Pm{z) 5(z - L)). 

The mean height h of the brush is the first moment of the probability density P{z)', the MC 
resuhs for the ratio h/Rg are plotted in Fig. Q vs oj^"^^''^^ for L = 200, 400 and 800. The data 
show a clear cross-over between two linear regimes with different slopes around Og ~ 10. The 
hnearity agrees with the prediction of Alexander-de Gennes scaling [|JJ, but the change in slope 
suggests a cross-over from the SCF regime to the regime dominated by strong excluded volume 
correlations which enhance stretching of the chains. Although the slope (~ 5.2) of the low Og 
regime is close to that of the earlier estimates obtained for different polymer models [l5l O, a 
clear-cut cross-over to the strong stretching regime has not been reported earlier. Indications in 
that direction are contained in the early work of Grest et al. [fT6l [TTl . By using a continuous 
model of polymer brushes, they observed deviations from the SCF scaling regime for grafting 
densities of 0.15 and chains 200 monomers long, resulting in enhanced stretching of the brush. 
However there was no clear evidence of a sharp transition between the SCF and the fully stretched 
scaling regimes, similar to to that shown in the end-monomer distribution function (Fig. [3]), or in 
the scaling plot (Fig. |4]) of the mean brush height. This was probably due to insufficient polymer 
length and to the presence of an attractive potential between the monomers. 

Let us know consider the "roughness" of the brush/solvent interface may be characterized by 

r -|l/2 

the relative fluctuation /m of the end-monomer altitude, (zij) — {zm) / {zm) which is readily 
derived from Pm{z)- The MC results for fM are plotted in Fig. |4] versus o for L=50 and 800. The 
difference is dramatic: while /m saturates rapidly at a value close to 0.4 for the shortest chains, the 
fluctuation is seen to decrease steadily for the longer chains, pointing to a relatively smooth upper 
interface of the brush. It is worthwhile to note that the corresponding fluctuation of the center of 
mass (CM) of the chains follows very similar curves, slightly below those for fj^- 

Sampling of stretched polymer conformations at high grafting densities becomes increasingly 
difficult for long polymers, such that Og > 10^. A way out is to switch to a multi-blob represen- 
tation of the grafted chains, as recently proposed for homopolymers and block copolymers in the 
bulk lfT8l[T9l . Within this coarse-graining procedure each chain is divided into n blobs of length 
I = L/n and radius of gyration Vgr^V. The average blob density within the volume of the brush of 
height is pft = Nn/{Ah) ~ A^L/z^^^^^^/^^, while the overlap density of blobs is ~ 1 /r^ ~ 1 //^v 
so that Pb/Pb ~ (L/na^/^^)^^^^. The minimum number of blobs required to ensure that blobs 
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Figure 2: (Color on-line)Same as in Fig. [TJ but for L = 200 and 0.002 < a < 0.25. Increasing deviations 
from SCF scaling regime are observed for a > 0.05 

do not, on average, overlap is hence n ~ o^^^. Effective interactions between the CM's of the 
bonded and non-bonded blobs and between a blob and the substrate, are determined by averag- 
ing over monomer degrees of freedom for given CM-CM distances. The coarse-grained multi- 
blob model requires effective pair potentials v(r) and (|)i (r) between the CM's of non-bonded and 
bonded blobs on the same or different (in the case of v(r)) chains, as well as an effective wall-blob 
potential and an effective tethering potential ^2{r) between the CM of the first blob and the 
grafting centre. All these effective interactions are determined by inverting MC results for the pair 
distribution functions between the CM's of a single grafted polymer made up of a small number 
of blobs [fTSl [T9l l20l . and are assumed to be transferable to finite grafting density conditions, as 
long as the system is in the weak overlap (p/, < pp regime. The form of the repulsive potential 
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Figure 3: (Color on-line)End-monomer height distribution function of grafted SAW polymers of length 
L = 800 versus z/L for grafting densities 0.002 < a < 0.25. Note the rather sharp transition to full stretching 
when a ~ 0.05, which contrasts with the gradual stretching observed for shorter chains in Fig. [T]j2] Inset: 
scaled monomer distribution function versus scaled altitude for the same values of a. Deviations from the 
SCF scaling regime set in for o > 0.01, and the profiles become nearly rectangular (full stretching) at the 
highest grafting densities. 



v(r) is Gaussian to a good approximation 

v(r) 
kRT 



■Ae 



-a{r/rgf 



(2) 



with A ~ 1.8 and a ~ 0.8. The tethering potential ^\{r) is well represented by the sum of v(r) 
and a harmonic spring potential, while ^2{r) is a similar harmonic spring potential, and the wall- 
blob potential is an exponential repulsion of range ~ Vg ||20|| . MC simulations of the multi- 
blob representation of grafted polymers sample single blob and total chain displacements with a 
standard Metropolis acceptance criterion, which are sufficient because of the softness of the above 
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Figure 4: (Color on-line)Plot of the mean height h =< zm > of a brush of SAW polymers (divided by 
Rg), versus G^g ^'^^^ for L = 50,200,400 and 800. The slopes of the two straight lines fitted to the low Og 
and high Og data indicate a cross-over from the SCF regime to the fully stretched regime around ~ 10. 

The slope (~ 5.2) of the low Og regime (dashed line) is close to earlier estimates obtained for different 

r 1 1/2 

polymer models Inset: relative fluctuation /m = {zj^) — {zm) / {zm) of the brush height versus 

v)/2v = 50 and L = 800, The fluctuations drops rapidly beyond the cross-over density, pointing to 
the homogeneity of the highly stretched brush. 

effective potentials. If n is chosen such as to satisfy the weak overlap condition above, the effective 
interactions can be extracted from low density MC simulations. This coarse-graining leads to an 
enhancement of sampling efficiency by several orders of magnitude, because of the reduction of the 
total number of degrees of freedom by a factor /, and the softness of the effective interactions. MC 
results for the CM probability densities, as determined for coarse-grained polymer chains of n = 5 
and 20 blobs are shown in Fig. [5} The profiles are seen to obey the SCF scaling at intermediate 
grafting densities, and to exhibit the enhanced stretching at very high grafting densities, in good 
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Figure 5: (Color on-line)Scaled monomer distribution function versus scaled altitude for brushes of multi- 
blob polymers of length n = 20 blobs, and 1.8 > a;, > 0.6 where Oh = or^. Increasing deviations from 
SCF scaling regime are observed for Ob > 1.2. Inset: comparison between the scaled monomer distribution 
function of a SAW polymer of length L = 800 at a grafting density o = 0.008 (triangles) and a multi-blob 
polymer under the corresponding physical condition of Oh = 0.617 (continuous curve). The plot shows very 
good agreement between the two representations. 

agreement with the full monomer- level MC results of Fig. [T]j4j 

In summary, we have shown by MC simulations of long grafted polymer chains that the brush 
profile and height undergo a rapid transition when the grafting density Og = oR^ > 10. Beyond 
that value, the probability density P{z) switches rapidly from universal, quasi-parabolic profile 
to a quasi-rectangular shape, while the fluctuations of the end-monomer altitude drops sharply, 
pointing to a rather flat brush/solvent interface. These changes are induced by correlations between 
monomers on neighbouring chains, which are neglected in SCF theory. We conjecture that these 
correlations will lead to an effective, entropic repulsion between stretched polymers, similar to the 
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Helfrich interaction [|2T1l between stacked membranes, which in turn may induce an ordering of the 
grafting centres on the substrate if these centres are mobile, as for copolymers anchored at a liquid- 
liquid interface, just as the Helfrich repulsion leads to a lamellar phase of stacked membranes. Our 
results point to new possibilities in experimental realization of chemically active surfaces, because 
the height distribution of the end monomer controls the chemical activity of the brush. Moreover 
the activity of such soft surfaces could be controlled by the length of the grafted polymers. 
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